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The influence of electron-electron Coulomb repulsion on nonsequential double ionization of rare- 
gas atoms is investigated. Several variants of the quantum-mechanical transition amplitude are eval- 
uated that differ by the form of the inelastic electron-ion rescattering and whether or not Coulomb 
repulsion between the two electrons in the final state is included. For high laser intensity, an entirely 
classical model is formulated that simulates the rescattering scenario. 



Multiple ionization of atoms by intense laser fields 
can proceed via different quantum-mechanical pathways. 
The atom or ion may be ionized step by step such that the 
transition amplitude is the product of the amplitudes for 
single ionization. If such a factorization is not possible, 
one speaks of nonsequential multiple ionization. Vari- 
ous physical mechanisms may be envisioned to underlie 
the latter, but electron-electron correlation is a necessary 
precondition. 

The actual presence of the nonsequential double ioniza- 
tion (NSDI) pathway was inferred long ago |lj from data 
at 1053 nm at rather low intensity, but the mechanism re- 
sponsible for it could only be identified after the reaction- 
microscope technique provided much more detailed infor- 
mation about the process than was available before 0- 
In principle, this technique is capable of analyzing double 
ionization in terms of all six momentum components of 
the ion and one electron, while earlier experiments were 
only able to yield total double-ionization rates. To the 
extent that the photon momentum can be neglected, this 
is synonymous with a complete kinematical characteriza- 
tion of the process. As a result, rescattering has emerged 
as the dominant mechanism, as it is for high-order har- 
monic generation and high-order above-threshold ioniza- 
tion |3J. For a review of recent developments, see Ref. |^|. 

An exact description of NSDI is, in practice, only pos- 
sible for helium and requires the solution of the six- 
dimensional time-dependent Schrodinger equation [3, Q . 
Various approximations and models have been consid- 
ered, in particular, one spatial dimension for each elec- 
tron 0, HI, or neglecting the backaction of the outer 
electron on the inner 0, or density- functional methods 
[ic| . Classical-trajectory calculations E2 yield reason- 
able agreement with the data far above the threshold. 
The exact quantum-mechanical transition amplitude has 
been analyzed in terms of Feynman diagrams [13. The 
basic diagram that contains rescattering has been eval- 
uated by several groups 

113, El EH It incorporates (i) 
tunneling of the first electron and (ii) its subsequent mo- 
tion in the laser field, (iii) the inelastic-rescattering col- 




FIG. 1: Feynman diagrams corresponding to the transition 
amplitude (a) without and (b) with electron-electron re- 
pulsion between the two electrons in the final state. The verti- 
cal wavy line and the dots in (b) indicate Coulomb interaction, 
which is exactly accounted for by the two-electron Volkov 
solution. The dashed vertical line represents the electron- 
electron interaction V12 by which the second electron is set 
free. 

lision with the second (up to this time bound) electron, 
and (iv) the propagation of the final two electrons in the 
laser field. It does not contain, for example, the inter- 
action of the electron in the intermediate state with the 
ion nor the interaction of the final electrons with the ion 
or with each other. 

In this Letter, we incorporate the Coulomb repulsion 
between the two electrons in the final state into the basic 
quantum-mechanical transition amplitude of NSDI, and 
compare with experiments. Moreover, we formulate a 
classical inelastic-rescattering model that reproduces the 
quantum-mechanical results for high laser intensity. 

The quantum-mechanical transition amplitude formal- 
izing the assumptions of the rescattering model is |l6j | 

/OO f-t 
dt / dt'(^ p2 (t)\V 12 
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xu[ v \t,t')v 1 \4 1 \t'))®\4 2 \t)). (i) 

It is pictorialized in the Feynman diagrams of Fig. ^ Ini- 
tially, both electrons are bound in their (uncorrelated) 

ground states |V4 (*')) = e^ Ean ^' IV'o )> where E an is 
the binding energy of the nth electron. At the time 
t', the first electron is released from the binding poten- 
tial Vi through tunneling ionization, whereas the second 



2 



electron remains bound. Subsequently, the first electron 
propagates in the continuum described by the Volkov 
time-evolution operator u[ V \t, t'), gaining energy from 
the field. At the later time t, it dislodges the second 
electron in an inelastic collision mediated by the interac- 
tion V±2 , which is accounted for in the lowest-order Born 
approximation. Throughout the paper, we compare two 
possible choices for this interaction: the Coulomb inter- 
action V12 ~ |i"i — r 2 | _1 and a three-body contact inter- 
action V12 ~ 5(ri — r 2 )<5(ri). The latter might be inter- 
preted as an effective electron-electron interaction on the 
background of the ion. For the final two-electron state 
l^pYpaW) with asymptotic momenta pi, p 2 , we take the 
correlated (outgoing) two-electron Volkov state 

l^(*)) = I^W)®l^(*)> 

x 1 F 1 (- l7 , l;t(|p||r| - p • r))e-^ 2 r(l + i 7 ), (2) 

which exactly accounts for their Coulomb repulsion 
[Fig. Gib)], and compare it with the product state of 
one-electron Volkov states \ippP(t)) [Fig. IHa)]. Here 
P = (Pi Pa)/2, r = n - r 2 , and 7 = 1/ (2|p|) (Coulomb 
repulsion is turned off by setting 7 = 0). 

In order to evaluate the multiple integrals in the tran- 
sition amplitude l|T|l. we expand the Volkov propaga- 
tor U[ v \t,t') = /d 3 k|^ V) (t))(Vk V) (*')l in terms of the 
Volkov states, (r\^\t)) = (2tt)- 3 / 2 exp{i[k + A(t)] • r} 
eocp[i5k(t)], with S k [t) = -(1/2) J* dr[k + A(t)] 2 the 
action of a free electron in the presence of the laser field 
described by the vector potential A(i). Since the Volkov 
solutions are plane waves (with time-dependent momen- 
tum), the spatial integrals yield the two form factors 

Vp k = ( P2 +A(0,Pi+A(t)|T/ 12 |k + A(t),4 2) ), (3) 
V k0 = (k + A(t')\V\4 1) ), (4) 

which, for the contact and Coulomb potentials that we 
consider, can be obtained in closed form. The remaining 
integrals over the intermediate-state momentum k, the 
ionization time t', and the rescattering time t are eval- 
uated by saddle-point integration, which is justified for 
the high ponderomotive energies Up of the experiments. 

The standard saddle-point approximation reduces the 
five-dimensional integration over k, t, and t' to a sum 
over the complex solutions k s , t s , t' s (s = 1, 2, . . .) of the 
saddle-point equations. The transition amplitude 

M (SPA) = -y (27rz) 5 / 2 Vp k3 14^ ^ gp(ta ^ ika) (5) 
s ^detS£(i,tf,k)| a 

is the coherent superposition of the contributions of all 
relevant saddle points. Here, the various exponentials 
in the amplitude have been collected into the ac- 
tion S p (t,t' The method is explained in detail in 



Refs. [2G3. In rescattering problems, the complex solu- 
tions come in pairs that approach each other very closely 
near the classical cutoffs. In this case, the standard 
saddle-point approximation becomes inapplicable. We 
here employ a so-called uniform approximation [Til Il9| , 
which invokes the same information on the saddles but 
works regardless of their separation. 

Both in the saddle-point and in the uniform approx- 
imation, the upshot of employing the correlated two- 
electron Volkov state is very simple: the result of 
the uncorrelated two-electron Volkov state just has to be 
augmented by the factor 27T7/[exp(27T7) — 1] [icj . 

The consequences are illustrated in the next figures, 
where Vyi is the contact interaction or the Coulomb in- 
teraction, respectively, and the parameters correspond 
to the data for argon (U H3 (note 0) and neon 0. 
We follow the presentation of the experimental data and 
decompose the final-state momenta into their compo- 
nents parallel and perpendicular to the (linearly polar- 
ized) laser field, so that p, = Cp<ii, P<_l) {i — 1,2). Then, 
we present density plots of the double-ionization proba- 
bility as a function of the parallel components pi\\ and 
P2\\, while the transverse components Pi± are partially 
or entirely integrated over. 

In general, final states where the two electrons 
have similar (vector) momenta are suppressed by their 
Coulomb repulsion, since the two electrons are set free 
simultaneously by the interaction Vn- In the (pi|i,f>2||)- 
correlation plots, this tends to reduce the population on 
the diagonal. The effect is most pronounced if both trans- 
verse momenta are small, and still sizable if only one is 
small [Figs. |21 El (c-f)]- When one transverse momentum 
is large, the effect is only moderate [Figs. El E] (a,b)]. 

Figures El El (a-d) are for the parameters of Ref. |2lj . 
While inconclusive if one of the transverse momenta is 
large [Figs. El El ( a ,b)], including the Coulomb repul- 
sion leads to improved agreement if one of them is small 
[Figs. El El (c,d)], but the resolution of the experiment 
does not allow one to settle either on Fig. Eld) or on 
Figs. H2c,d). Closer scrutiny of Figs. CJa)-(d) unveils a 
lack of symmetry with respect to the diagonal, which, 
possibly, has been observed in Ref. |2l|. It arises since 
the Coulomb form factor JSJ is only symmetric under in- 
terchange of all momentum components pi <-» p 2 [in the 
final state (0 , transverse and parallel momentum compo- 
nents can be interchanged independently]. Higher inten- 
sity and more restricted transverse momenta are investi- 
gated in Ref. , and for very small transverse momenta 
[Figs. El Et e ; f)] we do find the characteristic Coulomb 
pattern of Fig. El reflected in the experiment. The data 
agree better with Fig. EJe) , which does not incorporate 
final-state repulsion, than with Fig.EJf), which does. 

In neon |24j . the pjj_-integrated momentum correla- 
tion is quite well reproduced by the contact interaction 
|l6j . Corresponding results without and with final-state 
Coulomb repulsion are shown in Figs.EJa,b) andE^c,d), 
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FIG. 2: Comparison of the double-ionization probability den- 
sities without (left-hand column) and with (right-hand col- 
umn) electron-electron repulsion in the final state, as a func- 
tion of the electron momenta parallel to the laser field. The 
interaction V12 is specified by the three-body contact interac- 
tion. Parameters are for argon (E01 = 0.58 a.u., E02 = 1.0f5 
a.u.), the laser frequency is u) = 0.057 a.u. (Ti:Sa). Panels (a) 
and (b): I = 2.5 x 10 14 Wcm" 2 (U P = 0.54 a.u), |p ix | > 0.5 
a.u.; Ref. |2l|: (c) and (d): as before, but with |pu_| < 0.5 
a.u.; (e) and (f): I = 4.7 x 10 14 Wcm" 2 (Up = 1.0 a.u), |p 1± | 
or | P2 _l| < 0.1 a.u., Ref. I^. 



respectively. We find that taking into account the final- 
state Coulomb repulsion does not improve the agreement. 
The very simplest model - a contact interaction and no 
final-state repulsion - works best. This suggests that the 
V12 contact interaction can be interpreted as a reasonable 
zeroth-order effective interaction that takes into account 
the presence of the ion, which shields the long-range in- 
teraction between the returning and the bound electron. 

The rescattering diagrams of Fig. ^ though fully 
quantum-mechanical, stimulate a classical interpretation. 
In what follows, we recast their physical content into an 
entirely classical expression. First, let us enumerate the 
quantum features inherent in the transition amplitude 
||TJ: (i) the electron enters the field via tunneling, which 
enforces, in the saddle-point solution, complex t' s , t s , and 
k s ; (ii) the contributions of the individual saddle points 
are added coherently in the sum (in) in principle, the 
atom can absorb an arbitrary number of photons from the 
laser field, that is, double ionization still occurs, though 
at a much reduced rate, if the maximal kinetic energy of 



FIG. 3: Same as Fig. [21 but with V12 specified by the Coulomb 
interaction. 



the returning electron is less than I-E02I; (iv) the electron 
wave packet spreads during its propagation in the contin- 
uum. However, we do not expect the quantum features 
(ii) and (iii) to have a significant impact on the yields 
of double ionization far above the threshold, and feature 
(iv) is only relevant when long orbits contribute |25| . 

In order to check this surmise, we consider the classical 
model 



\M, 



class I 



dt'R{t')8 \E iet {t) - 
(P2 + A(t)) 2 



(Pi+A(i)) : 



2m 



2m 



E, 



02 



it; 



P k(t)i 



(6) 



Here, the electron appears in the continuum with zero 
velocity at the time t' at the time-dependent rate R(t'), 
which describes the quantum-mechanical tunneling pro- 
cess (cf. the above feature (i); we will use a simple tun- 
neling rate |2{J). For each ionization time t', the return 
times t = t(t'), and the corresponding kinetic energies 
E Tet (t) and drift momenta k(t) are calculated along the 
lines of the classical simple-man model 0. They corre- 
spond to the saddle-point solutions t s , except that they 
are real j^]. The <5-function in Eq. © expresses en- 
ergy conservation in the inelastic collision that sets free 
the second electron. The actual distribution of the fi- 
nal momenta is governed by the form factor |V r p k(t)| 2 - In 
Eq. the probabilities of the various orbits s = 0, 1, . . . 
are added so that their contributions cannot interfere, in 
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contrast to the quantum- mechanical amplitude (J3J). In 
Fig. 0{e) and (f), we present the result of the classi- 
cal model when the transverse momentum components 
Pi± are completely integrated over. They agree with the 
quantum- mechanical results of Figs.0{a) and^Jb) quite 
well. A similar conjecture was derived from the compari- 
son of one-dimensional classical-trajectory and quantum 
calculations jjj. However, the quantum distributions are 
slightly wider and the low-density regions are enhanced 
in comparison with the classical distributions, reflecting 
the fact that the quantum distributions extend into the 
classically forbidden region. 

In summary, we have investigated the effects of 
electron-electron repulsion in the final two-electron state 
of nonsequential double ionization. The calculations al- 
low us to conclude that footprints of electron-electron 
repulsion have been observed in experiments where the 
transverse momentum of one of the electrons is small, 
with better agreement when the repulsion occurs dur- 
ing the rescattering process rather than in the final state 
of the ionized electrons. We suggest that experiments 
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FIG. 4: Double-ionization probability densities for neon 
(E 01 = 0.9 a.u., E 02 = 1.51 a.u.), I = 10 15 Wcm" 2 (Up = 2.2 
a.u.), integrated over all transverse momenta |24|. Left pan- 
els: V12 the contact interaction, right panels: V12 the Coulomb 
interaction. The upper four are computed from the quantum- 
mechanical amplitude Q, in the absence (a,b) and presence 
(c,d) of final-state Coulomb repulsion, respectively. The re- 
maining panels (e,f) are calculated from its classical analog 
© without final-state electron-electron repulsion. 



in which both electrons are restricted to small transverse 
momenta would be most incisive. If the laser intensity 
is high enough, the quantum-mechanical momentum dis- 
tributions can be well reproduced in a purely classical 
model. 

Added note: Very recently, a paper by M. Wecken- 
brock et al., Phys. Rev. Lett. 91, 123004 (2003), was 
published, which addresses similar questions. 
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work was supported in part by the Deutsche Forschungs- 
gemeinschaft. 



* Present address: Quantum Optics Group, Institut fur 
theoretische Physik, Universitat Hannover, Appelstr. 2, 
30167 Hannover, Germany 
[1] A. PHuillier, L.A. Lompre, G. Mainfray, and C. Manus, 

Phys. Rev. A 27, 2503 (1983). 
[2] T. Weber et al., Phys. Rev. Lett. 84, 443 (2000); R. 

Moshammer et al, Phys. Rev. Lett. 84, 447 (2000). 
[3] PB. Corkum, Phys. Rev. Lett. 71, 1994 (1993). 
[4] R. Dorner et al, Adv. At., Mol., Opt. Phys. 48, 1 (2002). 
[5] D. Dundas, K.T. Taylor, J.S. Parker, and E. S. Smyth, 

J. Phys. B 32, L231 (1999). 
[6] H. G. Muller, Opt. Express 8, 425 (2001). 
[7] R. Panfili, S.L. Haan, and J.H. Eberly, Phys. Rev. Lett. 
89, 113001 (2002); S.L. Haan, P.S. Wheeler, R. Panfili, 
and J.H. Eberly, Phys. Rev. A 66, 061402(R) (2002). 
[8] M. Lein, E.K.U. Gross, and V. Engel, Phys. Rev. Lett. 

85, 4707 (2000). 
[9] J.B. Watson et al, Phys. Rev. Lett. 78, 1884 (1997). 
[10] D. Bauer and F. Ceccherini, Opt. Express 8, 377 (2001). 
[11] L.-B. Fu, J. Liu, and S.-G. Chen, Phys. Rev. A 65, 

021406(R) (2002). 
[12] A. Becker and F.H.M. Faisal, J. Phys. B 29, L197 (1996). 
[13] A. Becker and F.H.M. Faisal, Phys. Rev. Lett. 84, 3546 

(2000); ibid. 89, 193003 (2002). 
[14] R. Kopold, W. Becker, H. Rottke, and W. Sandner, Phys. 

Rev. Lett. 85, 3781 (2000). 
[15] S.V. Popruzhenko and S.P. Goreslavskii, J. Phys. B 34, 
L239 (2001). 

[16] S.P. Goreslavskii, S.V. Popruzhenko, R. Kopold, and W. 

Becker, Phys. Rev. A 64, 053402 (2001). 
[17] A. Becker and F.H.M. Faisal, Phys. Rev. A 50, 3256 

(1994). 

[18] C. Figueira de Morisson Faria and W. Becker, Laser 

Phys. 13, 1196 (2003). 
[19] C. Figueira de Morisson Faria, H. Schomerus, and W. 

Becker, Phys. Rev. A 66, 043413 (2002). 
[20] This is an excellent approximation and exactly true in so 

much as in the saddle-point integration the intermediate 

momentum k can be treated as real; for details, see C. 

Figueira de Morisson Faria et al. (to be published). 
[21] R. Moshammer et al, Phys. Rev. A 65, 035401 (2002). 
[22] M. Weckenbrock et al, J. Phys. B 34, L449 (2001). 
[23] In argon, as opposed to neon, there is experimental ev- 



5 



idence of an additional ionization mechanism, possibly 
excitation of the second electron followed by tunneling 
[B. Feuerstein et al, Phys. Rev. Lett. 87, 043003 (2001)], 
which generates momenta around pu\ — p 2 || = 0. These 
are absent in our calculations, which only include instan- 
taneous inelastic rescattering. The different behavior of 
neon and argon has been traced back to the respective 
ionization and excitation cross sections by V.L. Bastos 
de Jesus et al. (in print). 

[24] R. Moshammer et al., J. Phys. B 36, L113 (2003). 

[25] Spreading dampens the contributions of the longer or- 



bits, whose influence on the yields is minor, except at 
channel closings; see S.V. Popruzhenko, P.A. Korneev, 
S.P. Goreslavski, and W. Becker, Phys. Rev. Lett. 89, 
023001 (2002). Only the two shortest orbits are included 
in our calculations. 

[26] L.D. Landau and E. M. Lifshitz, Quantum Mechanics 
(Pergamon, Oxford, 1977). 

[27] In fact, the classical simple-man model for NSDI has been 
derived from the amplitude Q in Ref. [Tflj . 



